アレ イ ・ ア ン テ ナ を 用 いて 電波 の 到来 方 向 
が 分 か る 


到来 方 向 推定 シス テム の 中 人 と 実装 例 


Minseok Kim 


これ まで アレ イ ・ ア ン テ ナ を 用 いて 所 望 の 信号 を うま く 合 成す を 悪化 させ る と し て 問題 と な っ て いま 図 1). 
る 方 法 に つい て 説明 し た . し か し , 実は これ 以外 に も さま ざま この た め , 不法 リ ピー タ を 検出 する 装置 と し て , 複数 の 
な アプ リケーション が ある . 本 章 で は , アレ イ ・ ア ン テ ナ を 用 地点 に お いて 不法 電波 の 到来 方 向 DOA : Direction of 
いた 電波 の 到来 方 向 推定 法 (Direction of Arrival Estimat Arrival) を 測定 し , その 位置 が 特定 で きる シス テム が す で 
ion) に つい て 解説 し , FPGA を 用 いた 回 路 実装 を 紹介 する . に 開発 され て いま す . また , アレ イ ・ アンテナ を 用 いて 推 

(筆者 ) 定 する 手法 も いろ いろ と 検討 され て いま す . 
携帯 電話 に よる 緊 急 時 の 通報 が 増加 し て いま す . 米国 連 
_ 電波 の 到来 方 向 推定 0 FCC 0 緊 
(DOA Estimation) 生 M8H0 特定 で きる 機能 Enhanced 
91{ E-911) の 制度 化 を 行い まし た . 上 
最近 , 携帯 電話 が 目覚 まし い 勢い で 普及 し て いま す . し が 望ま れ て いま ずり リ . や は りこ の 場合 に も , アレ イ ・ アン 
か し , 建物 の 中 な ど に は 電波 の 届か な い 場 所 が 多く 存在 し , テ ナ を 人 
対策 と し て リピ ビー 人 多 Repeater) と 呼ば れる 電波 中 継 器 を 利 識 さ れ て いま ず 図 2). 

用 し て いる 場合 が 多く あり ます . ここ で , 誰 可 を 受け て い さら に , 図 3 の よう に セル ラ 移 動 通信 に お いて , ユー ザ 
な い 違 法 な リピ ビー タ が 使用 され た 場合 , これ ら の 放射 する 端末 の 到来 方 向 情報 を 有効 に 利用 する こと が 考え られ ます . 
不法 電波 が ノイ ズ と な り , 携帯 電話 や その ほか の 通信 環境 例え ば , 所 望 ユ ー ザ の 方 向 に は 送信 電力 を 集中 させ , ほか 
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図 1 不法 中 継 器 の 問題 図 2 救急 シス テム に お ける アレ イ ・ ア ン テ ナ の 応用 
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の ユー ザ 方 向 に は 指向 性 の ヌル を 形成 する こと で 干渉 を 低 
減 で きま す . 搬送 波 周波 数 が 上 下 回 線 で 異な る 周波 数 分 割 
複 信 FDD : Frequency Division Duplex ) に お いて 有効 
な 方 法 と し て 知ら れ て いま す . 到来 方 向 情報 は 上 下 回 線 に 
お いて 統計 的 に 等 し いこ と が 分 か っ て いる の で , 上 り 回 線 
( 移動 端末 っ 基地 局 ) で 推定 し た 到来 方 向 情報 を 下り 回 線 
( 基地 局 一 移動 端末 ) に 用 いる こと が で きま す . ここ で は , 
伝搬 路 の 変化 に 追い つけ る よう に 高速 で ,。 か つ , 高 精度 の 
到来 方 向 推定 処理 が 必要 に な り ま す . 

この よう に, 電波 を 用 いた 無線 局 端末 の 位置 推定 は , そ 
の ほか に も レー ダ ・ システム , 伝搬 路 同 定 シ ステ ム か ら 
近年 で は さま ざま な 無線 通信 シス テム や セン サ ・ ネ ッ ト 
ワー ク , コグ ニテ ィ ブ 無線 まで , その 応用 が 広がり つつ あ 
る と いえ ます . 

アレ イ ・ ア ン テ ナ を 用 いた 到来 方 向 推 定 手法 は , フー リ 
工 変換 の 原理 に 基づい て いま す . 最も 基本 的 な 手法 の ビ 
ム ・ フ ォ ー マ 湾 Beamformer), 線形 予測 渡 Linear 
Prediction) に 基づく 手法 と し て Capon 法 や 最大 エン ト ロ 
ビ 法 , さら に , アレ イ 入 力 信号 の 相関 行列 の 固有 値 展 開 に 
基づく 手法 と し て 多重 信号 分 離 決 以下 MUSIC: Multiple 
Signal Classification) や ESPRIT( Estimation of Signal 
Parameters via Rotational Invariance Techni ques) 法 
が 考案 され まし た . また , マル チ パ ス の 影響 を 受け や すい 
環境 に は チャ ネル 特性 の デー タベース 化 に よる 指紋 法 
( Finger Print) な ど が 提案 され て いま す . 

これ に つい て も っ と 詳細 な 理論 が 知り た い 方 は , 参考 文 
献 2 を 参考 に し て くだ さい . 

本 章 で は , アレ イ ・ ア ン テ ナ に お ける 到来 方 向 推定 法 に 
つい て , その 基礎 原理 を 説明 し , 現在 , 最も 幅広 く 利用 さ 


2 
- で る 7 


の 


反射 ・ 散 乱 波 較 
( a) 移動 通信 に お ける 伝搬 状況 凶 
図 3 セル ラ 移 動 通信 に お ける 到来 方 向 推定 の 利用 


ワイ ヤレ ス 通 信 の 剖 率 を 高め る 


誠 叶 


れ て いる 多重 信号 分 離 淫 MUSIC) に つい て 解説 し ます . 
また , FPGA を 用 いた ハー ド ウェ ア 構 成 と 信号 処理 部 の 実 
装 例 を 簡単 に 紹介 し ます . 


2. 到来 方 向 を 推定 する 方 法 


到来 方 向 推定 手法 に つい て まとも に 説明 し よう と する と , 
非常 に 専門 的 な 話 に な っ て し まい ます . ここ で は , その 原 
理 を 直感 的 に 捉え られ る よう に , で きる 限り 簡単 に お さら 
い を し ます . 

アレ イ ・ ア ン テ ナ に よる 到来 波 の 方 向 推定 に は , 古い 歴 
史 が あり ます . 指向 性 ビー ム に よる 空間 走査 を 行う ビー 
ム ・ フ ォ ー マ 法 が 最も 簡単 な 手法 と 知ら れ て いま す . ビー 
ム ・ フ ォ ー マ 法 は , 一 様 励 振 アレ イ ・ アンテナ の メイ ン ・ 
ロー ズ 主 ビ ー ム ) を 全 方 向 に 対し て 走査 し , 得 ら れ た 電力 
スペ クト ル か ら 出 力 電力 が 大 きく な る 方 向 を 探索 する 方法 
で ず 2. 

図 4 の よう に , 重み 係数 を ある 角度 9 に する た め に は , 
第 2 章 の 説明 より 以下 の よう に アレ イ ・ ス テア リン グ ・ ベ 
クト ル に 設定 する こと で , メイ ン ・ ビ ー ム を 角度 の に 向け 
る こと が で きま す . 


eo=|1 mp っ 学 ene| ーー 
ん 

HERE 《 1) 

ee っ 学 ew-Dmd|| 


ここ で , 4 は アン テ ナ 素 子 数 , 2 は アン テ ナ 素 子 間隔 を 
示し ます . この と き , アレ イ 出 力 は 以下 の よう に な り ま す . 


y( の =@T『(6) + FFP ドキ ( 2②) 


的 と する 信号 較 


に ささ 


( b) 到来 方 向 か ら 所 望 ユ ー ザ に 電力 を 集中 させ る こと が で きる 図 
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共役 


ここ で , 記号 *※ は 共役 を 示し ます . 

ビー ム ・ フ ォ ー マ 法 は , 空間 領域 で の フー リ 工 変 換 で 関 
係 付け られ , FEFT( Fast Fourier Transform) を 用 いて 高 
速 処理 が 可能 に な り ま す . この よう な 電力 検出 方 法 で は , 
理想 的 な が ビーム は ある 角度 に 対す る イン パル ス 的 な 応答 を 
持ち ます . 現実 的 に は 制限 され た アレ イ ・ ア ン テ ナ 素子 の 
数 に より , メイ ン ・ ロ ー ズ メイ ン ・ ビ ー ム ) の 半値 幅 が 広 
く な り ま す . また , サイ ド ・ ロ ー ブ か ら も 信号 を 受信 し て し 
まう ため, 一般に 角度 分 解 能 が 悪い こと が 知ら れ て いま す . 
し か し , シス テム 的 に は 安定 に 動作 し , 最も 簡単 に 実現 で き 
る こと か ら , 現在 も 幅広 く 用 いら れ て いま ず 2. 

ディ ジタル 制御 が 自由 に で きる よう に な っ て か ら , さま 
ざま な 手法 が 考案 され て きま し た . その 中 で , MUSIC や 


Sv 入射 波 較 


( a) リニア ・ ア レイ ・ ブ ロッ ク 図 図 


- 40Li H i i H i 
-80 -60-40-20 0 20 40 60 80 
角度 ) 
( b) 8 素子 アレ イ ・ ア ン テ ナ に +30' 方 向 か ら 入 射 さ れる と き の 図 
電力 スペクトル 図 


図 4 ビー ム ・ フ ォ ー マ 法 の 原理 


114 Design Wave Maggzine 2007 December 


ESPRIT と いう 部 分 空間 法 を 用 いる 超 高 分 解 能 手法 が 考案 
され , 現在 , 最も 注目 され て いま ず 9 ?. 


3. 多重 信号 分 離 法 


(MUSIC アル ゴリ ズム : MUltiple 
Slgnal Classification) 


江 衣 の よう に た ビーム: ウォー ヤマ 油 で は メイツ ヽ ビーム 
の 半値 幅 に よっ て 角度 分 解 能 が 決ま る の で , 到来 波 の 細か 
い 分 離 は で きま せん . MUSIC アル ゴリ ズム は , アレ イ 入 
カ デ ー タ か ら 作成 し た 自己 相関 行列 の 固有 値 展開 に より , 
信号 部 分 空間 と 雑音 部 分 空間 に 分 離 で きま す . それ ら が お 
互い に 直交 する と いう 面白 い 性 質 を 用 いる こと で , 直 高 分 
解 能 の 到来 信号 の 推定 精度 が 実現 で きま す . 本 章 で は 少し 
専門 的 な 話 に な り ま す が , その 原理 に つい て 科 略 に 説明 し 
ます 。 

MUSIC アル ゴリ ズム に お ける 信号 処理 は 図 5 の よう に な 
り ま す . アレ イ ・ アンテナ に 入力 さ れ た 信号 は , 受信 器 で 
周波 数 変換 と AD 変換 され た 複素 ディ ジタル 信号 に な り ま 
す . 各 素 子 で の 雑音 を 含ん だ 受信 信号 ベク トル /) か ら 以 
下 の よ うに 自己 相関 行列 を 作成 し ます . 


し 2 ト て OSI 議 Pom 《 ? 


ここ で , 琴 ・] は 期待 値 を 示し ます . この 相関 行列 は , 
信号 成分 と ノイ ズ 成 分 と し て , 以下 の よう に 分 離し て 考 
えら れ ま す . 信号 ベク トル & の の 相関 行列 を Ass, 複素 
白色 ガウ ス ・ ノ イズ を M ) と する と , ノイ ズ 相 関 行 列 は 
その 分 散 値 と 単位 行列 か ら o 47 と な り , xx は 以下 の よ 
う に 変形 で きま す . 


0 お 00 ま : の 3 すっ 0 
= ん ル cs Le ター ーー ーー 《 5) 


実際 に 信号 と ノイ ズ の 分 離 は 固有 値 分 解 に より 行わ れ ま 
す . 一 般 に 相関 行 多 xx は 非 不定 値 エ ルミ ー ト 行列 
( Positive Definite Hermitian Matrix ) に な り , 固 有 値 分 解 
に より ベク トル 空間 を 次 式 の よう に 信号 部 分 空間 と ノイ ズ 
部 分 空間 と し て 分 離 で きま す . 


Az =0。 個 償 ro 人 物 時 ーー…… (6) 


ここ で の ひ 。 と ひ , は, それぞれ 信号 と ノイ ズ の 部 分 空間 
の 固有 ベク ト ル で す . 4 。 は 対 角 成 分 が 実数 の 固有 値 { ん 


1 ん ル 2, "っ ん 衣 ん 生ん 2 生 … 三 えん) の 対 角 行 列 で , 
4 y は 対 角 成 分 を o 2 と する 対 角 行 列 に な り ま す . ノイ ズ 
部 分 空間 の 固有 ベク ト ル は , 信号 部 分 空間 の 固有 ベク ト 
ル と 直交 する 性 質 が あり ます . 
全 方 向 に 対す る ステ アリ ング ・ ベ クト ル と ノイ ズ 固 有 べ 
クト ル の 内 積 を 求め ます . この 値 が 角度 スペ クト ル に な り 
ます . 直交 する 二 つ の ベク ト ル の 内 積 は 0 と な る た め , ス 
ペク ト ル の 値 が 0 に な る 角度 を 電波 の 到来 角度 と みな せま 
す . これ は , 鋭い ヌル を 利用 する こと に な り , 高 分 解 能 の 
E 定 が 可能 に な る わけ で す . これ が メイ ン ・ ビ ー ム の 操作 
に よる ビー ム ・ フ ォ ー マ 法 と は 大 きく 違う 点 で す . 

計算 し た スペ クト ル の スプ リ アス を 抑え る た め に , 一 般 
こ は 複 数 の ノイ ズ 固 有 ベ クト ル を 用 いて 同じ 処理 を 行い , 
スペ クト ル を 平均 化し ます . また , 計算 され た スペ クト ル 
値 の 逆数 を と る こと で , 図 5 の よう に ヌル 点 が 負い ピー ク 
と し て 現れ る こと に な り ま す . 図 5 の スペ クト ル は , 4 素 
子 ア レイ ・ ア ン テ ナ を 用 いて 三 つ の 信号 の 到来 方 向 を 推定 
し た 結果 と な り ま す . 


5 ge"( の ・ (の ) 

7 g7(@ の ). 役 Uz(⑥) 

=|er 1…ek|: ノ イズ 固有 ベク トル 凶 の 
( 素子 アレ イ に 信号 
が / 波 到来 する 場合) 


4. 到来 方 向 推定 シス テム 


ここ で , 横浜 国立 大 学 で 試作 され た 到来 方 向 推定 シス テ 
ム に つい て 紹介 し まず 5? の. 推定 手法 は , 複素 数 行列 演 
算 を 実数 の み の 演 算 に 低減 で きる 行列 変換 ユニ タリ 変換 
と いう ) を 用 いた ユニ タリ MUSIC 法 で す . MUSIC 法 は , 
前 述 の よう に 非常 に シン プル な 原理 で す . 単純 繰り 返し 処 

に 適し て いる た め , FPGA を 用 いた 高速 演算 が 可能 と 考 
えら れ ま す . 


人 @ シス テム の 構成 

ディ ジタル 信号 処理 回 路 は , 大 規模 FPGA 米国 Altera 
社 の Stratit EP1S40」) に 実装 し まし た . シス テム の 構成 
を 図 6 に 示し ます . ハー ド ウェ ア は 第 2 章 で 説明 し た も の 
と 同じ で す . 到来 方 向 推定 に お ける 信号 処理 は , 相関 行列 
の 計算 , 固有 値 分 解 処理 , FFT に よる スペ クト ル 計 算 , ス 


ペク ト ル 探 索 と いっ た , 大 きく 四 つ の 処理 プロ セス に 分 け 
られ ま す . 

まず , 8 素子 り ニア 
RF 信号 5GHz) を , 受信 器 で 周波 数 変換 し , AD コン バー 
タ で 直接 サン プリ ング 可能 な 低い IF 周波数 40MHz) に し 


・ ア レイ ・ アンテ ナ か ら 受 信 さ れ た 


ます . AD コン バー タ に より サン プリ ング し , 12 ビ ッ ト の 
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スペ クト ル を 計算 図 


MUSIC スペ クト ル 凶 
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5 MUSIC 法 の 概念 的 な 計算 フロ ー 
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に 


ー タ と し て FPGA に 入力 し ます . この 信号 


ディ ジタル ・ 
2 変換 され , 複素 
数 の 解析 信号 に な り ま す . 

検 濾 後 の 解析 信号 は 一 般 に オー バ サ ン プリ ング され て い 
る た め , サン プル 数 を 減ら す 必 要 が あり ます . 各 素 子 か ら 


の 信号 を ベク ト ル と し て 取り 扱い , 自己 相関 行列 を 作成 し 
ます . ここ で , 内 部 処理 は 16 ビ ッ ト の 固定 小数 点 演算 で 行 
われ ます . それ か ら , 相関 行列 の 固有 値 展開 を 行い , 雑音 
固有 ベク ト ル を 計算 し ます . 
固有 値 展 開 は 周期 ヤマ コ ビ 淫 Cyclic Jacobi) を 用 いま し だ 9. 
これ は ベク ト ル 回 転 の 単純 な 繰り 返し の み で 実装 で きま す . 
ベク ト ル 回 転 動作 を 並列 化す る こと が 容易 な た め , FPGA 
上 の 回 路 化 に も っ と も 適し て いる よう に 思わ れ ま す . 角度 
スペ クト ル の 計算 は , 空間 領域 で の FFT を 用 いる こと に 
より 高速 処理 が 可能 で す . その スペ クト ル 値 に 対し て 1 次 
元 の ピー ク 検索 を 行う こと で , 到来 角 を 推定 する と いう 構 
成 に な り ま す . 

信号 処理 は AD 変換 基板 上 の FPGA で 行い ます . シス テ 
ム 制御 や 外部 イン ター フェ ー ス は , SH4 CPU 上 で NetBSD 


受 
信 
機 
図 


が 動作 する 組み 込み CPU 基板 を 使用 し まし た ( 9. 


@ 信号 処理 の FPGA 実装 
ここ で , 今回 設計 し た 信号 処理 部 の FPGA へ の 実装 に つ 
いて 簡略 に 説明 し ます . 誌面 の 関係 で 詳細 な 説明 は 省略 し 
て し ます が , 興味 の ある 方 は 参考 文献 を 参照 し て くだ さい . 
ユニ タリ MUSIC 法 に お ける 信号 処理 の 流れ を 表 1 に 示 


し ます . 第 一 ステ ッ プ は , ユニ タリ 変換 で 入力 信号 を 実数 
化す る こ NN RE IE IE 
行い ます . 8 素子 の アレ イ の 場合 に 実数 の 相関 行列 は 
7 2 78 
= 紹 6 SH RDG ( 8) 
781 782 “| 788 


の よう に な り , 各 要 素 は 次 の 式 で 計算 され ます . 


エタ br 


や 


刀 | (な 7 の 7 りー o の 7o ⑦) 
の Yo の 7/( の キャ 7(2x/o⑰) 


引 拉 り 


に 
に 
= 


図 6 

FPGA 上 に Unitary MUSIC 
アル ゴリ ズム を 実装 し た シ 
ステ ム の ブロ ッ ク 図 


相関 行列 了 尼 16 
計算 較 | 
M 
スペ クト ル 狗 
の 図 
ピー ク 検 出 図 
8 
SH4 図 
Ethemnet | トー ペー 
: CPU ボー ド 図 


処理 内 容 図 


ベク ト ル に 対し て ユニ タリ 変換 を 行う 


空間 平均 処理 と 時 間 平 均 処 理 を 行う 較 


CORDIC を 用 いる Cyclic Jacobi 固 有 値 分 解 を 行う 凶 


演算 式 較 
ユニ タリ MUSIC に お ける 主 ニー の 
な 演算 防 
&、⑦⑰ = 記 (w-D+(- め が. 例 Rely の の )| 
国 
屯 = と V の I&,| 
さ 2 
>|pr7te 血 | 


と +1 


ノイ ズ 固 有 ベ クト ル の ユニ タリ 逆 変 換 , 図 
FFT に よる 角度 スペ クト ル 計 算 図 
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リス ト 1 
対象 性 を 用 いて 対 角 


に VHDL に よる 


記述 を 
の 28 要 素 の 乗算 回 路 を 並列 で 行う こと に な り ま す . FPGA 


示し ます . 64 要 素 の 中 , 


の 専用 乗算 ブロ ッ ク を 用 いる こと で , 高速 な 回 路 実 装 が 可 
能 に な り ま す . 作成 され た 相関 行列 を 用 いて , 固有 値 展開 


リス ト 1 8 素子 の アレ イ 入力 信号 か ら 相関 行列 を 計算 する 回 路 corrxx.vhd) 


ーー 素 天天 天天 天天 天天 天天 天天 天天 広末 素 天天 天天 天天 天天 天天 大 天天 天天 天天 ーー 


Compute Corre1at1on Mat エ ュ 1x 
Corrxx .YhQd 


ーー 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 天天 大 天天 天天 大 ーー 


ーー Rxx Rea1 ( 
-- 16 Snapshotg 
ーー 12b1E Tnput8g, 


Y*Herm1 モ ian (Y) 
16bitg Outputg 


LTBRARY TEEE: 

USE TEEE.std 1ogio 1164 .a11: 
USE TEEE.std 1og1io ARTTH.a11: 
USE TEEE.std 1og1o STGNED .a11 
USE work.u1] paokage .ALh : 


ENTTTY Co エエ xx TS 
PORT ( 
で 1k 
で 1kx 
ena 
rege 上 


TN 
TN 
TN 
TN 


Std 1og1o: 
std 1og1o: 
gtd 1og1o: 
std _1og1o: 


TN 
TN 


81Y12a 了 : 
81Y12a 了 : 


-- 12bitg 4 array 


OUT sg1Y16ar36) : 
END ENTTTY CO エエ XX : 


ARCHTTECTURE エ て ヒ 上 ] OF Corrxxx エ 8 


COMPONENT un1Earyr TS 
PORT ( 
で 1k 
ena 
rege 


エ N std 1og1o: 
エ N std 1og1o: 
エ N std 1ogio: 


T エ N 81V12a エ 了 : 
TN 81v12a エ 了 : 


ーー 12 bitg 3 array 


OUT 81Y12a エ : 
Yd OUT s1V12a 了 ) : 
END COMPONENT un1 ヒ a エ ャ ヒエ : 


ーー 12 bitg 3 arraY 


COMPONENT boxcar TS 
PORT ( 
で 1k 
rege 上 
ena 


エ N std 1og1C: 
エ N std 1og1C: 
エ N std 1ogio: 
x_in エ N std 1ogio_ vector(15 downto 0): 
Y_ou OUT std 1og1c vector(15 downto 0) ) : 
END COMPONENT boxoar 了 : 


COMPONENT 11rupdate TS 
PORT ( 
で 1k 
rege 
ena 


TN 
TN 


std _1og1C: 
Std 1og1G: 
エ N std 1ogio: 
x 1n エ N std 1ogio veoctor(15 downto 0): 
Y_ou OUT std 1og1o vector(15 downto 0) ) : 
END COMPONENT 11rupdae: 


STGNAL 
STGNAL 
STGNAL 
STGNAL 
STGNAL 
STGNA ェ L 
BEGTN 


11 , 
YYy11, 
yy1 
窟 


tg1 
yyq1 


81Y12a 了 : 
81Y12a エ 了 : 
81V12a エ 36 
81Y16a エ 36 
81Y16a エ 36 
Std 1og1C: 


in 
mu1t x, mu1t y 
8 


ェ , reg_ エ 


を 行い ます . 

ハー ド ウェ ア 実 装 の 容易 な Cyclic Jacobi 法 を 用 いま す が , 
回 路 は CORDIC 法 に より 固定 小数 点 演算 が 実現 で きま ず 5?. 
固有 値 分 解 に より 得 ら れ た ノイ ズ 固 有 ベ クト ル を 複素 数 に 


戻し て , その 離散 フー リ 変換 に より 角 


度 ス ペク トル を 計 


-- INPUTs 
ti1 (1 TO 8) 
tq1 (1 TO 8) 


<= 1 xi1(1 TO 8) : 
<= 1 xq(1 TO 8) : 


un1 モ a エ ャ ヒエ 1 ユ 
tg1, yy11, 


umn1 ヒ ary キ エエ PORT MAP 
YYyq1) 


(clk, '1! , rege 七 , 


も 1 。 


1oop_r ]]: 
2 人 
Fixradd (Fixxmu1t (yyi1 (]) , yy11 (] ) 
Fixmu1t ( 


FOR ] TN 1 TO 8 GENERATE 
) <= 
,16) , 
yyq1 (}) , yyq1 (]) , 16) ) : -- r]]} 
END GENERATE: 
1oop_r 1]: FOR ] TN 2 TO 8 GENERATE 
エエ (]+ 7) <= 
Fixadd (Fixmu1t (yyi1 (1) , yy11 (]) , 16) , 
Eixmu1Et (yyq1 (1) , yyq1 (] ) , 16) ) : -- エ 1] 
END GENERATE: 
1oop_r 2]: FOR ] TN 3 TO 8 GENERATE 
エエ (]+13 ) <= 
fixadd (Fixmu1t (yyi1 (2) , yy11 (]) , 16) , 
Eixmu1t (yyq1 (2) , yyq1 (]) , 16) ) : -- エ 2] 
END GENERATE: 
1oop_r 3]: FOR ] TN 4 TO 8 GENERATE 
エエ (]+18 ) <= 
fixadd (F1xmu1E (yy11 (3) , yyi1 (]) , 16) , 


Fixmu1t (yyq1 (3) 


,YYd1 (] ) , 16) ) , -- エ 3] 
END GENERATE: 
1oop_r 4]: FOR 〕 TN 5 TO 8 GENERATE 
エエ (]+22 ) <= 
Fixadd (Fixmu1t (yyi1(4) , yy11 (]) , 16) , 
fixmu1t (yyg1 (4) , yyq1 (}) ,16) ) : -- エ _4] 
END GENERATE: 
1oop_r 5]: FOR 〕 TN 6 TO 8 GENERATE 
エエ (]+25 ) <= 
Eixadd (Fixmu1t (yyi1 (5) , yy11 (]) , 16) , 
fixmu1t (yyg1 (5) , yyq1 (}) ,16) ) : -- エ _5] 
END GENERATE: 
1oop_r 6]: FOR 〕 TN 7 TO 8 GENERATE 
エエ (]+27) <= 
ixadd (Fixmu1t (yyi1 (6) , yy11 (]) , 16) , 
fixmu1t (yyq1 (6) , yyg1 ( け ) ,16)): -- F 6] 
END GENERATE: 
エ (36) 
fixadd (Fixmu1t (yyi1 (7) , yy11 (8) ,16) , 


Eixmu1t (yyq1 (7) , yyq1(8) ,16) ) : 


<= 
2) 


1oop0 : FOR 〕 TN 1 TO 36 GENERATE 
PROCESS (cl1k, ena, reset) T8 
BEGTN 
エ F re8e 七 THERN 
in エ (}) (OTHERS=>!0') : 
ELSTF olk'event AND c]1k = !1! THEN 
in エ (]) <= エエ (]) 』 
END TTF: 
END PROCES8: 
END GENERATE: 


「]「 


<= 


-- end 1oop0 


1pf 1oop0 : FOR 1 TN 1 TO 36 GENERATE 
1pfFuu 11rupdate PORT MAP 
(c1k, reset, ena, 1n て (1) , 
END GENERATE 1pfF ]oop0: 


ェ ( せ ) ) 』 


END ARCHTTECTURE エ ヒ ]: 


Design Wave Maggzine 2007 December 117 


表 2 ユニ タリ MUSIC の FPGA 
実装 結果 準 同期 検波 処理 を 含む ) 
アン テ ナ 素 子 数 
4 素子 
( 空間 平均 処理 付き ) 
4 素子 
8 素子 


MM 
前 依 / 
送信 アン テ ナ 財 " 


| 


LM 人 


( a) 電波 暗室 に お ける 実験 の 様子 較 
ユニ タリ MUSIC 法 で 計算 し た 角度 スペ クト ル の 実践 


図 7 


算 し ます . フー リエ 変換 は FPGA ベン ダ 提 供 の IP 
( Intellectual Property) コ ア な ど を 用 いれ ば 手軽 に 実装 で 
きま す . この 角度 スペ クト ル か ら ピー ク 点 を 検索 する こと 
で , 到来 方 向 が 推定 され ます . 

表 2 に Altera 社 の StratiX EP1S40) を 用 いて 実装 し た 結 
果 を 示し ます . 並列 処理 を する こと を 考慮 し 実装 面積 は や 
や 大 きめ と な り ま し た が , 8 チャ ネル の 準 同期 検波 を 含め 
て 8 素子 アレ イ ・ アンテ ナ の 実装 まで 1 チッ プ に 収まり ま 
し た . 

7 a) は 吸収 体 に より 反射 が 抑え られ る 電波 暗室 に お け 
る 実験 の 様子 示し ます . 8 素子 アレ イ ・ ア ン テ ナ に - 5" 
と 20' か ら そ れ ぞ れ 2 波 が 到来 する 場合 の 結果 を 図 メ b) に 
示し ます . 到来 方 向 が 正しく 推定 され る こと が 分 か り ま す . 


へ 


の へー へ 


と 


へ 


ンス 


人 @ 終わ り に 

今回 の 特集 は , 最近 の 移動 通信 の 動向 や MIMO など で 注 
目 を 集め て いる アレ イ ・ ア ン テ ナ 信号 処理 技術 の 基礎 原理 
に つい て 専門 外 の 人 に も 優し く 解説 し , FPGA を 用 いた 実 
装 例 も 紹介 し まし た . 少し 専門 的 な 話 に な っ て し まっ た と 
ころ も あり まし た が , 興味 の ある 方 は ぜひ , 参考 文献 を 読 
ん で 頂き た いと 思い ます . 

これ まで 1 対 1 で 扱わ れ た 情報 が , 複数 の スト リー ム で 
伝送 され , 信号 処理 に より 合成 され る よう に な り ま し た . 
より 高速 か つ 並 列 な 処理 が 必要 に な っ て くる 中 , FPGA 特 
有 の 信号 処理 能力 は と て も 魅力 的 で ある と 思い ます . 無線 
通信 分 野 に お ける FPGA の 活躍 を 期待 し た いと 思い ます . 


T 
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角度 ”) 図 


( b) 8 素子 アレ イ ・ ア ン テ ナ に 2 波 が 図 
それ ぞ れ - 5 と 20" か ら 到来 する 図 
場合 の 実験 結果 較 
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